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We consider spin- half quantum antiferromagnets in two spatial dimensions in the quantum limit, 
where the spins are in a valence bond solid (VBS) phase. The transitions between two such VBS 
phases is studied. In some cases, an interesting second order transition controlled by a fixed line 
with varying critical exponents is found. A specific example is provided by an antiferromagnetically 
coupled bilayer system on the honeycomb lattice where a continuous quantum phase transition can 
generically exist between two VBS phases. Furthermore, these critical points are deconfined, in the 
sense that gapped spin-f /2 spinon excitations emerge right at the transition. The low energy physics 
of this critical point (upto marginally irrelevant interactions) contains just a free quadratically 
dispersing 'photon'. The phase structure on one side of this continuous transition is very intricate 
consisting of a series of infinitely closely spaced further transitions in a 'devil's staircase' form. 
Analogies with previous examples of deconfined quantum criticality are emphasized. Closely related 
transitions in single layer systems are explored. These are second order only at some multicritical 
points. The solvable Rokshar-Kivelson point of quantum dimer models of single layer systems is 
found to correspond to a non-generic multicritical point 



I. INTRODUCTION 

Although the theoretical study of quantum phase 
transitions^ has recieved much attention over the last 
several decades, it remains a source of rich and unex- 
pected physics to this day. For example, recent work 2,3 
on the quantum criticality of certain quantum magnets 
(spin-1/2 quantum antiferromagnets on the square lat- 
tice) has shown that a direct and generically continuous 
transition between the Neel state and a valence bond 
solid state is possible. Such a transition violates 'Lan- 
dau's rules' of classical phase transitions which prohibit a 
direct continuous transition between phases that posess 
such different symmetries. Perhaps more interestingly 
the critical theory is unusual and naturally expressed in 
terms of new emergent, 'deconfined' degrees of freedom 
that carry fractional spin, along with an emergent U(l) 
gauge field. This is despite the absence (i.e confinement) 
of such fractional spin excitations or the associated gauge 
field in either phase. Such critical points were dubbed 
'deconfined' quantum critical points (QCPs). In this pa- 
per we shall examine another example of a "deconfined" 
quantum critical point which occurs in a different context 
i.e. between two VBS phases. 

The existence of such deconfined QCPs may be surpris- 
ing since the gauge theories arising in condensed matter 
systems (apart from physical electromagnetism) are gen- 
erally compact, i.e. are defined in terms of a periodic 
gauge connection e M rather than a single valued non- 
compact vector potential a. In D=2+l, such gauge the- 
ories are well known to be "always" (but see below) con- 
fining. In particular, the simplest and apparently generic 
compact U(l) gauge theory hamiltonian is 
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with integer valued "electric field" E and 27r-periodic vec- 
tor potential a (canonically conjugate to one another) de- 
fined on the links of a two-dimensional lattice. Here A is 
the lattice gradient, and Ax is the lattice curl, defined 
as the gauge flux through a plaquette. In the natural 
cases arising in condensed matter contexts, one imposes 
the Gauss law constraint A • E = (the "even" gauge 
theory) or A • E r = e r on bipartite lattices, with e r = ±1 
taking opposite signs on the two sublattices (the "odd" 
gauge theory). In both cases, regardless of the values of 7 
or the "dielectric constant" < e < 00, the gauge theory 
is always confining. The value of the nominally "irrele- 
vant" (in the renormalization group sense) coupling K is 
then immaterial. 

Remarkably, Ref. |2j concluded that a deconfined phase 
is possible if a gauge theory similar to above is coupled to 
gapless matter fields, which arise naturally in quantum 
antiferromagnets at particular QCPs. In this paper we 
study a different class of transitions which provide other 
instances of 'deconfined' quantum criticality in quantum 
magnets. Although in contrast to Ref. 0, they are not 
examples of Landau forbidden transitions, they neverthe- 
less display several interesting similarities to those priors. 
We consider transitions between two valence bond solid 
(VBS) states: quantum paramagnets with a gap to spin 
excitations. In either phase these spin-carrying excita- 
tions are conventional (though gapped) spin-1 magnons 
or their composites. We study direct second order tran- 
sitions between two such phases where 'deconfinement' 
obtains right at the critical point. We show that the 
critical mode can be viewed as a gapless deconfined U(l) 
gauge field with a quadratic dispersion. Moreover the 
magnon excitations, that are sharply defined gapped ex- 
citations on either side of the transition, break up into 
two weakly interacting (gapped) spin-1/2 spinons. These 
spinons are minimally coupled to the critical deconfined 
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U(l) gauge field. This results in a weak interaction be- 
tween the spinons so that they are essentially free. 

In contrast to the examples of deconfined quantum 
criticality of Ref. which are strongly interacting^, the 
critical fixed points discussed here have a simple free field 
description. In particular, they may be considered as a 
limit of the compact gauge theory in Eq. JQ| where the 
dielectric constant e — > oo. At this point, the energy cost 
of an emergent electric flux to quadratic order in E is pro- 
portional to | A x E| 2 rather than |E| 2 as usual. Despite 
the absence of gapless matter fields, this can provide a re- 
alization of a deconfined QCP, provided the system has 
sufficient (physically achievable) symmetries, as shown 
below. We emphasize that the appropriate physical sym- 
metries occur naturally, for instance in an antiferromag- 
netically coupled spin-1/2 honeycomb bilayer. Moreover, 
the QCP in this case is "generic" , i.e. only one physical 
parameter needs to be tuned (corresponding to tuning 
1/e through zero) to reach it. While particular "quantum 
dimer" models (see below) have provided some concrete 
realizations of compact U(l) gauge theories, many other 
microscopic mechanisms leading to Eq. ^ exist, and it is 
fruitful to regard the gauge theory as the most general 
progenitor of such deconfined QCPs. 

Closely related phenomena have been mentioned pre- 
viously in the literature 5 ^ in the context of work on the 
quantum dimer model 7 . The dimer model provides a 
caricature of spin-gapped quantum paramagnetic phases 
and admits on certain lattices a solvable point known as 
the Rokhsar-Kivelson (RK) point. On a square or honey- 
comb lattice the RK point is known to be critical and to 
separate two phases with very different VBS order. An 
important development was the field theoretic descrip- 
tion of this RK point, that was conjectured by Henley 5 , 
and further elaborated in 6 . In this paper we will analyze 
the generic behaviour expected in such a transition be- 
tween VBS phases, and show that it differs in important 
ways from the physics of the RK point. In particular for 
single layer quantum spin systems on the square lattice 
we argue that generically there is no second order transi- 
tion between two VBS phases. It has been argued in Ref. 
that a similar situation also obtains for the single-layer 
honeycomb lattice, namely there is generically no second 
order transition. In contrast as we show below such a sec- 
ond order transition with 'deconfined' criticality obtains 
on the bilayer honeycomb lattice with antifcrromagnetic 
interlayer coupling. In fact, a line of fixed points with 
continuously varying exponents is obtained. The univer- 
sal long distance physics at one of these critical points 
corresponds to the long distance physics at the RK point 
(upto logarithmic corrections arising from a marginally 
irrelevant operator). Thus, although in general there are 
significant differences from the generic transition, the RK 
point in some cases can still provide useful information 
about the universal critical properties of these generic 
transitions. A unique feature of the RK point is that it 
occupies a very special place in the phase diagram of the 
generic spin model, a point that will be further discussed 



below. 

We now summarize our results for the continuous 
quantum phase transition between two VBS phases of 
spin-1/2 quantum antiferromagnets on the bilayer hon- 
eycomb lattice. The precise nature of the neighbour- 
ing phases themselves will be described shortly. There 
is a simple Gaussian description of the critical theory 
which is parametrized by a 'stiffness' K and has a dy- 
namical scaling exponent z = 2. In addition there is 
a marginally irrelevant operator that will lead to loga- 
rithmic corrections. Apart from the usual 'thermal' op- 
erator needed to tune to the critical point, the theory 
is found to have no other relevant perturbations for a 
range of stiffnesses, leading to a fixed line with contin- 
uously varying exponents depending upon K. Indeed, 
this is an interesting example of a fixed line in D = 2 + 1 
dimensions, and is closely related to three dimensional 
statistical physics models of the Lifshitz point in certain 
liquid crystal systems^. 

Our results rely crucially on a 'dual' description of the 
quantum paramagnetic phases of antiferromagnets (with 
collinear spin correlations) in terms of a sine-Gordon (or 
'height') field x^2*ii. It is constructed (see Section [Bj 
by defining Ei — EijAjX, so that all physical properties 
are invariant under the global shift x ~~ * X + 1 • This for- 
mulation may be obtained in a number of different ways 
which are briefly discussed in Section [B] below. In the 
continuum limit appropriate near critical points between 
various paramagnetic phases the Euclidean action of this 
model reads (for the honeycomb lattice bilayer system): 

S = So + Si + Si ns t (2) 

5 = \ Jd 2 xdr{(d T x) 2 + p(V X ) 2 + A'(V 2 x) 2 }(3) 

51 = Jd 2 xdr- A \A X \ i + ... (4) 

Sinst = - Jd 2 xdr Acos(27rx) (5) 

This continuum action is invariant under global integer 
shifts of the x field, as well as \ ~ * — X (E — > — E) 
which arises from symmetry under bilayer exchange. The 
ellipses represent other terms (higher derivatives and 
higher powers) which could be added that are consistent 
with the symmetries above. As discussed below, the field 
X may be interpreted as the dual of a U(l) gauge field. In 
this interpretation the A term describes instanton events. 
The transition of interest occurs when p changes sign 
with p = at the critical point. Instantons are relevant 
in either phase but will turn out to be irrelevant at the 
critical point for some range of K . In addition the quartic 
term w(Vx) 4 will be shown to be marginally irrelevant. 
Thus, the action Sq in Eqn. [21 (which was recently studied 
along with several interesting generalizations in Ref>i£) 
will describe, upto marginally irrelevant terms, the fixed 
point action of a generic quantum critical point between 
two VBS phases of the bilayer honeycomb lattice antifer- 
romagnet. When p > the sine Gordon field (also known 
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as the height field) x has zero 'tilt' in the ground state. 
It describes a featureless paramagnetic state of the origi- 
nal quantum magnet, which may be caricatured as being 
made of singlets formed on the interlayer rungs. How- 
ever when p < one expects a finite 'tilt' i.e. Xr ~ Q • r - 
Thus the transition may be expected to occur between a 
VBS phase with zero tilt of the height field and a VBS 
phase with a nonzero tilt of the height field - the tilt in- 
creasing 'continuously' from zero on moving away from 
the 'deconfined' critical point. Indeed while this general 
expectation is more or less correct, the detailed picture 
is somewhat more complicated owing to the relevance 
of the instanton operators away from the transition. In 
fact, as described in section IlII CI interesting structure 
obtains on the 'tilted' side of the transition. The back- 
ground electric field (tilt) changes from its value of zero 
at the critical point to a non-zero value some distance 
away, through a fractal sequence of interleaving regions 
of constant tilt (confined phases) and non-constant in- 
commensurate regions (deconfined phases) of zero width 
but finite total measure. The resulting structure is known 
as an 'incomplete devil's staircase' £ 

In contrast to the above picture, the RK point sits be- 
tween a VBS phase with zero tilt of the height field, and 
the staggered phase which has the maximum tilt to the 
height field, i.e. with 2irQ above one of the six recipro- 
cal lattice vectors of minimum length. These phases are 
illustrated in Fig. 2] One can then ask whether the RK 
point occupies a special place in the phase diagram of 
the generic model. In fact, the well known property of 
the RK point that the ground state in each tilt (winding 
number) sector are all degenerate indicates that indeed 
it occupies a special place in the phase diagram. We 
have already noted that the RK point exhibits the same 
universal critical properties as the generic model after 
tuning a few relevant and marginal operators. However, 
to reach the special position occupied by the RK point 
in the phase diagram, fine tuning of several dangerously 
irrelevant operators is required. These operators do not 
affect critical properties but determine the phase struc- 
ture in the vicinity of the critical point. The precise po- 
sition of the RK point in the phase diagram is dicussed 
below. Note, the generic phase structure can be recov- 
ered if we add to the RK point an operator that (for 
example) corresponds to the quartic term S± in equation 
(J2J). Moreover, this particular operator will also gener- 
ate logarithmic corrections to the RK point correlation 
functions, which will then precisely match the one of the 
generic critical points for the bilayer honeycomb system. 
Thus, generic behaviour, both in terms of critical prop- 
erties as well as phase structure neighbouring the critical 
point can be obtained by adding this one operator to the 
RK point Hamiltonian. 

In the usual quantum dimer model, this fine tuning of 
the Hamiltonian to access the RK point is achieved by 
keeping only terms that involve a single plaquette. How- 
ever, if we regard the dimer model as an approximate de- 
scription of some underlying quantum spin model, then 



in general we need to include dimer kinetic and potential 
terms on arbitrarily large-sized loops (but with coeffi- 
cients that decrease with increasing loop size). The ex- 
act degeneracy mentioned above then obtains only if the 
dimer kinetic energy on every such loop is set equal to the 
corresponding potential energy. Thus for an underlying 
quantum spin model this exact degeneracy presumably 
requires 'infinite' fine-tuning. 

A different way to draw a sharp distinction between the 
RK point and a generic (two parameter tuned) multicrit- 
ical point that allows for a direct zero tilt to staggered 
(maximum tilt) state transition, is to ask - what phases 
can be accessed from these points by a small change in 
bare parameters? From the generic multicritical point 
(shown in Figure |5Jd) three states - the staggered, zero 
tilt and an infinitesimally tilted state - can be accessed. 
From the RK point however, in addition to the above 
three states, states with arbitrary values of the tilt can 
be accessed (Figure |5J;) . This is a consequence of the 
exact degeneracy of states with different tilts at the RK 
point which gives rise to the special phase structure in 
its vicinity. 

The layout of this paper is as follows. In Section [H] we 
discuss the sine Gordon representation of this problem 
from different points of view, the non-linear sigma model 
approach of Haldane, the easy plane deformation of the 
spin half magnet of Lannert et and the familiar 

height representation of the quantum dimer model. In 
Section 11111 we analyze the spin half quantum antifcrro- 
magnet on the bilayer honeycomb lattice using this rep- 
resentation, and find a stable fixed line controlling the 
transitions. Various properties of these critical points 
are discussed, and the 'devil's staircase' phase structure 
is obtained. In section [V] the single layer square and 
honeycomb lattice spin-half antiferromagnet is discussed, 
where the generic transition between VBS phases is first 
order. In section IVII we discuss the well known RK 
points, and how they fit into the general structure de- 
scribed in this paper. 



II. SINE GORDON REPRESENTATION FOR 
PARAMAGNETIC PHASES OF COLLINEAR 
QUANTUM ANTIFERROMAGNETS 

Our analysis relies crucially on a formulationifi*ii of 
the physics of quantum paramagnetic phases in terms of 
a sine Gordon field theory on the dual lattice. In this 
Section, we describe this sine Gordon description and 
discuss its origin and interpretation from several different 
perspectives which together provide considerable insight. 
We will first discuss the single layer case and then move 
onto the double layer. The lattice Euclidean action for 
the sine Gordon model appropriate to a single layer has 
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the following structure: 

S = Si + S 2 + (6) 

Sl = \fY,^rXr) 2 + P {^X) 2 +K{^ X f (7) 

/oo 
dr A w cos (27rn(x r - a r )) (8) 

r n=0 

Here r runs over the sites of the dual lattice. The sym- 
bol A refers to a lattice derivative. The ellipses repre- 
sent other terms that are consistent with all the sym- 
metries (both internal and lattice) that could be added. 
In the second term n is an integer that runs from 
to oo . The a r are independent of time but vary in a 
definite manner on spatial lattice sites. On the square 
lattice, a r = 0,1/4,1/2,3/4 on four sublattices. Thus 
exp (2«7ra r ) oscillates rapidly on four sublattices. For a 
spin model on the honeycomb lattice, the correspond- 
ing sine Gordon theory is defined on the dual triangu- 
lar lattice. In this case a r = 0,1/3,2/3 on the three 
sublattices of the triangular lattice so that exp (2iira r ) 
oscillates rapidly on three sublattices. 

Near the phase transitions of interest in this paper, the 
universal physics is adequately described by a continuum 
limit of this action. The oscillating phase factors due 
to the a r means that in the continuum the only terms 
that survive from S2 are those with n = 0(mod4) for 
the square and n = 0(mod3) for the honeycomb lattices 
respectively. 

The physical basis of the sine Gordon description may 
be understood in many ways. We will mainly discuss the 
square lattice - the honeycomb lattice is very similar. 

Gauge theoeretic/quantum dimer description: 
First consider a description of VBS phases in terms of 
quantum dimer models. A dimer is to be considered a 
caricature of a single valence bond and is taken to live 
on the bonds of the original lattice. For spin-1/2 sys- 
tems, it is natural to constrain the dimer Hilbert space 
by requiring that there is exactly one dimer emanating 
out of each lattice site. It has long been appreciated that 
such quantum dimer models on bipartite lattices can be 
fruitfully viewed as compact U(l) gauge theories. This is 
understood as follows: First divide the bipartite lattice 
under consideration into A and B sublattices. To each 
dimer we may associate an integer-valued vector 'electric 
field' E that starts from an A sublattice point and ends 
at a B sublattice point. In terms of these electric fields, 
the dimer constraint simply becomes the Gauss law 

A ■ E = ±1 (9) 

where the + sign is for the A sublattice and the — sign 
for the B sublattice. Thus we see that the dimer Hilbert 
space is identical to that of a particular compact U{\) 
gauge theory with fixed background 'charges' ±1 on the 
two sublattices. Such theories were christened 'odd gauge 
theories' in Ref. |(j 



Alternately such a gauge theoretic description of the 
disordered phases of the antiferromagnet can be directly 
obtained by starting with a slave particle {eg. Schwinger 
boson or fermion) description. Consider for instance a 
Schwinger boson representation. This has a U{1) gauge 
redundancy associated with arbitrary phase rotations of 
the bosons at each site. In a spin gapped paramagnetic 
phase the Schwinger bosons can be integrated out and 
the physics is described by a compact U{1) gauge theory. 
In this route too an odd gauge theory obtains. 

There is a well-known duality mapping between com- 
pact U(l) gauge theories and sine-Gordon field theories 
- see for instance Ref. The interpretation of the sine 
Gordon field is as follows. In the absence of the com- 
pactness in the gauge theory, the total magnetic flux 
is exactly conserved. There is a corresponding (topo- 
logical) global U(l) symmetry. In the sine-Gordon rep- 
resentation this becomes an ordinary global U(l) sym- 
metry that corresponds to an arbitrary global shift of 
the sine-Gordon field. This symmetry is present in the 
sine Gordon model if the coefficient A„ of all the cosine 
terms is set to zero. Including compactness in the original 
gauge theory allows for 'instanton' events which destroy 
flux conservation. Indeed the flux can change in mul- 
tiples of 2tt. These are precisely captured by including 
terms like cos2-7m(xr — oc r ) in the sine Gordon descrip- 
tion. Thus e 27 ™* corresponds to an instanton event at 
which 2n units of gauge flux is created. The shift a r is 
due to the 'oddness' of the gauge theory. Specifically the 
presence of background charges in the gauge theory leads 
to Aharanov-Bohm phases for the gauge flux which are 
encapsulated in the shifts a r - for details see Re£i£. 

Semiclassical Description: Much further insight is 
obtained by a semiclassical perspective that is partic- 
ularly appropriate if there is considerable short-ranged 
Neel order. Deep in the Neel phase the long distance 
low energy fluctuations of the Neel vector are described 
by the familiar quantum 0(3) non-linear sigma model 
field theory. To correctly describe quantum paramag- 
netic phases it has been recognized for some time now 
that this continuum field theory must be augmented by 
appropriate Berry phase terms that are sensitive to the 
microscopic spin at each lattice site and to the details 
of the lattice structure. The Berry phase terms vanish 
for all smooth configurations of the Neel field. In two 
spatial dimensions note that such smooth configurations 
allow for topologically non-trivial configurations known 
as skyrmions. However as shown by Haldane the Berry 
phases are non-zero in the presence of singular configura- 
tions -known as hedgehogs or monopoles - in space-time. 
At the location of the monopole the skyrmion number as- 
sociated with the Neel field configuration changes. That 
such skyrmion tunneling events arc allowed is a conse- 
quence of the presence of a lattice in the microscopic spin 
model. The calculation of Ref. Ill shows that the Berry 
phases associated with a single monopole (defined on the 
plaquettes of the original lattice) oscillates from one pla- 
quette to another with amplitude exp 2ina r where a r is 
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as defined above on the dual lattice sites (or equivalently 
on the plaquettes of the original lattice). 

The quantum paramagnetic state is associated with a 
proliferation of such monopole events. The Berry phases 
associated with a single monopole event leads to it trans- 
forming non-trivially under lattice symmetry operations. 
Thus proliferation of single-strength monopoles leads to 
broken lattice symmetry in the paramagnetic state. 

To describe the different possible paramagnetic phases, 
it is convenient to imagine integrating out all gapped 
spin-carrying excitations and focus on the physics of 
the skyrmion fields and the associated monopole events. 
First consider a limit in which the monopoles are ignored 
(i.e imagine tuning the monopole fugacity to zero). In 
this limit the skyrmion number is exactly conserved. This 
corresponds to a hidden (topological) global U(l) sym- 
metry in the absence of monopoles. The paramagnetic 
phase may then be thought of as a condensate of these 
skyrmions so that this global U{\) symmetry is sponta- 
neously broken. The low energy excitations are fluctua- 
tions of the phase of the skyrmion field and will be gap- 
less. Indeed the corresponding action may be identified 
with Si above with x/27r being the phase of the skyrmion 
field. Including monopole events leads to explicit break- 
ing of this global £7(1) symmetry. Clearly in this picture 
the S2 term corresponds precisely to skyrmion creation 
events (which are the monopoles) with the appropriate 
Haldane oscillating phases encapsulated in the shift fields 
a r . 

Easy Plane Limit: Finally we briefly mention a 
derivation of the action above in the easy plane limit 
of the original spin model. As argued in Refiiiii (see 
Ref j£ for a physical discussion) , a very useful continuum 
description of easy plane spin- 1/2 magnets is provided 
by focusing on vortex fields (known as 'merons') in the 
(XY-like) order parameter. The dual action takes the 
form 

C = Lii, + C v + La + C-inst (10) 
= E \ {d» -iA ti )iP a \ 2 + r\rP\ 2 +u{\^\ 2 )\ll) 

o=l,2 

Cv = vlViH^I 2 (12) 

Ca = K(e^ K dvA K ) 2 (13) 

dust = -A 4 (W>i> 2 ) 4 + c.c) . (14) 

Here ipi% represent the two meron fields that are mini- 
mally coupled ( as is usual) to a non-compact U(l) gauge 
field A, and \ip\ 2 = IV'il 2 + |^2| 2 - As described in Ref.0 
for the square lattice the last term physically describes 
the monopole or instanton events discussed above for the 
isotropic models. This continuum model has a global 
Z4 symmetry associated with ±(7r/4 + mir/2) shifts of 
the phase of V>i and tp2 respectively (m = 0, 1,2,3). In 
this description valence bond solid phases correspond to 
(tpi) = (ip%) 7^ 0. In such a condensate, the gauge field 
A acquires a mass by the usual Anderson-Higgs mecha- 
nism and may be ignored at low energies. Furthermore 



the global Z4 symmetry is also broken - thus a low en- 
ergy description is provided by focusing on the relative 
phase 6 between tpi and ip2- Clearly the theory has the 
same structure as the continuum limit of Eqn. E|and we 
identify \ = 2ixB. This discussion readily generalizes to 
the honeycomb lattice - the main difference is that the 
monopole events are tripled (leading to Z3 symmetry). 

^From any one of these perspectives it is clear that 
(anti)vortices in the \ field correspond to spin-1/2 spinon 
configurations in the original spin model. Specifically we 
define a (anti)vortex as a point in space around which 

J dl- Vx = ±l (15) 

For instance in the dimer model these correspond pre- 
cisely to points where the dimer constraint is violated 
(i.e to monomers). Equivalently we note that skyrmions 
and spinons see each other— as sources of 2ir flux - so 
that a spinon configuration corresponds to a vortex in 
the skyrmion phase. 

Bilayer Systems: The discussion above is readily 
adapted to bilayer systems. Specifically consider a bilayer 
spin-1/2 quantum antiferromagnet on a square or honey- 
comb lattice. The symmetries of the microscopic Hamil- 
tonian now include the Ising-like layer exchange symme- 
try in addition to SU (2) spin rotation, time reversal and 
all the lattice symmetries. This layer exchange symmetry 
will play an important role in our analysis. Consider first 
the limit in which the interlayer antiferromagnetic ex- 
change on each rung is the largest coupling. In this limit 
it is appropriate to first diagonalize the 'rung' Hamilto- 
nian. For each rung, the ground state is a singlet and the 
first excited state is a triplet. A useful model of such a 
bilayer is to replace each rung by an 0(3) quantum rotor 
with the Hamiltonian 

H =\ T,^- J E ni-hj + (16) 

i <ij> 

Here hi is a unit three component vector defined on each 
rung (labelled by i) and Li is the corresponding angular 
momentum. The ellipses represent other short-ranged 
terms consistent with the symmetries. The rotor vector 
hi and the angular momentum Li have the same sym- 
metry properties as the difference and sum of the two 
microscopic spins on the rung at i respectively. Under 
layer exchange we then have 

hi — > —fii (17) 
Li -» Li (18) 

Thus layer exchange symmetry implies that the rotor 
Hamiltonian be invariant under the full group 0(3) of 
rotations (which includes improper rotations). In addi- 
tion time reversal is a separate symmetry that is imple- 
mented by an antinunitary operator that changes the sign 
of both h and L. We remark that this must be contrasted 
with single layer Heisenberg spin magnets. These can 
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also be modelled as quantum rotors but with appropri- 
ate monopoles placed at the origin of the n-sphere at each 
site. Now the presence of the monopoles implies that the 
improper rotations of the rotor vector are no longer sym- 
metries. Thus these must be regarded as 5*0(3) rotors. 
This distinction will also be important for us below. 

Consider now paramagnetic phases of the bilayer 
model. We will specifically be interested in phases that 
obtain close to the strong interlayer exchange limit where 
the 0(3) rotor description becomes appropriate. As with 
the single layer systems discussed above it will be con- 
venient to obtain a gauge theoretic description of these 
paramagnetic phases. This may be obtained by passing 
to a CP 1 description of the rotors in terms of spinon vari- 
ables z. The z fields are minimally coupled to a compact 
U(\) gauge field a but unlike the single layer case the 
mean spinon number is zero per site. In a mean field de- 
scription of paramagnetic phases the spinon fields will be 
gapped. Beyond mean field, integrating out the gapped 
spinons leads to a compact £7(1) gauge theory. The ulti- 
mate fate of the spinons is determined by whether or not 
this gauge theory is confined. Again in contrast to the 
single layer case the Gauss law constraint of this gauge 
theory is simply 

A-E = (19) 

with no background charges. Here is the 'electric' 
field defined on the links of the honeycomb or square 
lattice. As usual this is conjugate to the gauge field a^: 

K- ,£«]=». (20) 

What are the symmetries of this gauge theory? Clearly 
all the symmetries of the square or honeycomb lattice 
that forms each layer are also symmetries of the gauge 
theory. In addition the symmetry of the rotors under 
improper rotations (the layer exchange symmetry) im- 
plies that the gauge theory must be invariant under the 
discrete symmetry 

Eij -> -Eij (21) 
^ a>ij (^^) 

This may be seen in several ways. For instance we note 
that the magnetic field corresponding to the gauge field 
a is precisely the skyrmion density associated with the 
configuration of the n fields. The latter is odd under im- 
proper rotations of n (for instance there is a well-known 
expression for the skyrmion density as a trilinear in n). 
Similarly the electric fields correspond to the skyrmion 
currents which are likewise odd under improper rotations. 
Parenthetically we note that under (the antiunitary) time 
reversal the electric field is even while the gauge field is 
odd. 

This compact U(l) gauge theory is readily dualized to 
obtain a dual sine-Gordon description in terms of the \ 



field. The dual action takes the form 

S = Si + S 2 + (23) 

& = ^/ T E( 9 *) 2 +PN 2 +^(A 2 X) 2 (24) 

/oc 
dr YJ YJ A„ cos (27TOXr) (25) 

r n— 1 

There are two important differences with the single layer 
case. First the absence of background charges in the 
gauge theory implies that there are no offsets a r for the 
height fields Xr- Second the discrete layer exchange sym- 
metry implies that the action must be invariant under 

Xr > Xr- 

It is interesting to contrast the bilayer with a spin- 
3/2 antiferromagnet in a single honeycomb layer. In the 
latter the gauge theory appropriate to the paramagnetic 
phase may be viewed as an 'even' gauge theory, i.e one 
where there are no background charges in the Gauss law 
constraint. But nevertheless as the microscopic model 
is not invariant under improper rotations of the spin (or 
equivalently the rotor vector in a rotor description) the 
gauge theory does not have the discrete symmetry of Eqn. 
PHI associated with changing the signs of both E and a. 

III. THE BILAYER HONEYCOMB LATTICE 

In this Section, we specialize to the bilayer honeycomb 
lattice, assuming the presence of strong interlayer antifer- 
romagnetic coupling. We will consider the lattice valence 
bond solid phases and phase transitions of this spin half 
quantum antiferromagnet using the sine Gordon descrip- 
tion Eqn. 1231 It will be convenient to explicitly write out 
the lowest order non-linear terms that are allowed by the 
symmetries. These take the form 

S int = J dT^\A x \ 4 + ^Re[(A x + iA y ) x ] 6 (26) 

so that the full action is 

S = Si + S 2 + S mt (27) 

with S12 given in Eqn. We have included in Si nt 
the v term which gives the lowest order effect of the six- 
fold broken rotational symmetry of the hexagonal lattice. 
Note that it appears only at sixth order in Ax, and hence 
is nominally irrelevant at the critical point. It is never- 
theless important for p < (see below). Note that the 
symmetries of the bilayer system include X r —> ~Xr at 
each site r which corresponds to layer exchange. This 
symmetry forbids the appearance of terms that are odd 
in x- Such terms are allowed in single layer systems that 
do not have this symmetry - indeed Ref. identifies a 
cubic term that drives the transition first order. 

Let us for a moment consider how the various VBS 
phases arise from the action (|27|l . For p > 0, the system 
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(a) 



(b) 



FIG. 1: Caricature of VBS phases on the bilayer honeycomb 
lattice, (a) The zero tilt state, with singlet bonds (thick lines) 
on the interlayer rungs. Note, this state does not break any 
lattice symmetry, (b) One of six possible maximally tilted 
(staggered) phases. 



would like to have zero tilt A\ on the average. Moreover, 
we can ask what the effect of the monopole tunneling 
term is, the first such term is the one that inserts a single 
monopoles (n = 1) of Eanl25l and it is easily seen that 
this operator has long ranged correlations if p ^ 0, and 
is hence a relevant perturbation. The resulting phase 
will be one where the height field is pinned at a uniform 
value, and this state may be caricatured as one where the 
spins form singlets with their partners in the other layer 
as shown in Fig. ^l). 

For negative values of p, a state with a finite tilt is 
expected, i.e. 



k = (2ttAx> ^ 0. 



(28) 



Ignoring the monopole operators for the moment, the 
system will choose a tilt ko whose magnitude for small 
negative p is obtained at the mean-field level by minimiz- 
ing together the p and u terms, i.e. 



|kn|=27rW-^ 



(29) 



with corrections of 0[v /u{— p/u) 3 / 2 ] from v and other 
higher order terms. Fluctuation effects due to the 
marginality of u will slightly enhance |k | by a loga- 
rithmic factor of little importance. The direction of the 
vector ko is, however, determined by the sign of v. In 
particular, the six discrete directions with ko x + iko y = 
|k |e i2 ™/ 6 or k 0x + ik 0y = |k |e i277 ( m+1 / 2) / 6 , with m = 
. . . 5 are preferred for v < 0, v > 0, respectively. The 
added effects of monopoles, which will modify the true 
tilt vector to k =^ ko, however, will have to be carefully 
considered discussed in the last part of this section. 

In any case, a phase transition between a VBS with 
zero average tilt, and one with nonzero average tilt, cor- 



responds to taking p from positive to negative values. 
The phase transition that lies between requires that we 
look at the theory with p = 0. 



A. Stability of the Fixed Line Controlling the 
Transition 

The critical theory for such a VBS transition is then 
proposed to be: 

S <= f drJ2l{(drX) 2 + K(A 2 X f} (30) 

r 

we need to check that this simple Gaussian theory is 
stable against switching on a small monopole tunneling 
(term S% in the action l|23|0 and quartic interactions (Si n t 
of equation (|26|) ). Again, we consider correlators of the 
single monopole tunneling event, the n — 1 term. This 
now has the following power law decay in space: 



/ e i277 Xr (0) e -i27r X o(0)\ 



r Vk 



(31) 



which implies that the monopole insertion operator is 
irrelevant if > 4, or equivalently < K < (f) 2 - 

Thus, there is a line of fixed points, with different ex- 
ponents, that are parametrized by K and are all stable 
against switching on weak monopole tunneling. Indeed 
this is very similar to the line of fixed points obtained in 
D = 1 + 1 in a variety of systems such as the spin half 
XXZ chain. Note however that our theory has dynamical 
exponent z — 2. 

We now consider the effect of the quartic interaction 
term, Si n t in Eqn. 1261 which by naive power counting is 
marginal at these fixed points. Since we will be looking at 
values of K for which the monopole tunneling events are 
irrelevant, they are disregarded in the discussion below. 
We consider a continuum model of the critical theory 
(J20J, with modes restricted to wavevectors below a cer- 
tain cutoff A, We perform a one loop RG, assuming that 
we start with a small value of the interaction parame- 
ter u, and study its flow on integrating out the large 
wave- vector modes A(l — dl) < \k\ < A. After the appro- 
priate rescaling to keep the quadratic term (|3U[1 invari- 
ant, we obtain the following flow equation for the quartic 
coupling 9 : 



du _ 9 
dl ' 167rid 



(32) 



This implies that a quartic coupling with u > is 
marginally irrelevant, and the coupling flows back to zero 
logarithmically with distance. Therefore the long dis- 
tance physics in this case will be controlled by the critical 
action l|3U|l . with logarithmic corrections arising from this 
marginally irrelevant operator. Thus the critical points 
are stable towards turning on a quartic interaction for 
u > 0. (For u < however, the quartic coupling is rel- 
evant, and the transition is very likely driven first order 
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- as is indeed already the case in mean field theory for 
negative u). 



B. Properties of the Transition 

The analysis above has established the presence of a 
fixed line controlling the transition between VBS phases 
with zero tilt and those with a non-zero tilt of the height 
field. In this subsection, we discuss some properties of 
this fixed line. We first note that the irrelevance of the 
monopole tunneling terms implies that the global sym- 
metry of the continuum theory is enlarged to U(l). This 
corresponds to invariance under arbitrary global shifts of 
the height field \- I n terms of the U(l) gauge theory 
(whose dual is the sine Gordon theory) this implies that 
the compactness is asymptotically irrelevant all along this 
fixed line. Indeed as shown in Appendix 1X1 the free Gaus- 
sian action is readily seen to describe a quadratically dis- 
persing gapless photon in the gauge theory representa- 
tion. The gauge flux of the theory is then conserved. 
This signals deconfinement. However the monopoles are 
important to correctly describe the physics of the phases 
on either side of the transition - thus they represent 'dan- 
gerously irrelevant' perturbations. 

If we ignore the marginally irrelevant quartic term the 
fixed point theory has a free field form. Thus all as- 
pects of the critical behavior while non-trivial are emi- 
nently tractable. For instance it should be possible to 
compute real time dynamical correlators of, say, the op- 
erator e l2 ^ Xrl ^ at non-zero temperatures in the scaling 
limit. We will however not pursue this here. 

It is instructive to ask about the fate of the gapped 
spin-carrying excitations right at this critical fixed line. 
As mentioned in Section [HI the presence of a spin-1/2 
spinon at some spatial site leads to a vortex in the height 
field x- Away from the critical point, in (for instance) 
the zero tilt VBS, the relevance of the monople tunnel- 
ing terms lead to pinning of the height field which implies 
that there is a huge energy cost that increases linearly 
with system size for these vortices. More precisely con- 
sider a pair of spinons (on opposite sublattices of the orig- 
inal lattice). This generates a vortex-antivortex pair in \, 
The energy cost for separating this pair by a distance R 
grows linearly with R away from the critical point. Thus 
we have (linear) confinement of spinons and the elemen- 
tary spin-carrying excitations have spin-1. But right at 
the critical fixed line the enlargement of the symmetry to 
t/(l) implies that vortices in the x ne ld are cheap. Ignor- 
ing the quartic perturbation, an elementary computation 
shows that the energy cost of a vortex is finite indepen- 
dent of system size. Including the quartic term leads to 
a weak 1/R 2 interaction (up to logarithmic corrections 
due to the marginal irrelevance of u) between a vortex- 
antivortex pair separated by distance R - this then is the 
interaction between two spinons (on two opposite sublat- 
tices) separated by a distance R. Thus as expected the 
spinons are deconfined and free to propagate above a spin 



gap. 

Within the approximation of ignoring the weak interac- 
tion between the two spinons, the magnon spectral func- 
tion A(k,oj) at the gap edge is readily calculated. One 
finds a sharp step A(k, u) ~ 9(uj — A) where A is the spin 
gap. Thus the magnon spectral function has no quasipar- 
ticle peak and is anomalously broad. 

C. Devil's Staircase 

We now consider the behavior on the "tilted" side of 
the Lifshitz point, in which the compact gauge theory 
is expected to have a non-vanishing background electric 
field. The neighborhood of the Lifshitz point has been ar- 
gued in Ref y|to realize an "incomplete Devil's staircase" . 
While the essential features have already been sketched 
in Ref.Q, we recapitulate them here for completeness and 
to present a few additional points not mentioned therein. 
In particular, we will describe the thermal transitions of 
the commensurate tilted VBS states within the devil's 
staircase, and point out two distinct types of low energy 
excitations within these phases. 

Let us first think generally about the nature of phases 
with some finite and non-zero background electric field, 
i.e. not in the direct vicinity of the putative critical point. 
At this stage we do not take any continuum limit, work- 
ing with a field Xr defined on the discrete lattice sites (of 
the dual square or hexagonal lattices). If one neglects at 
first the compactness of the gauge theory, i.e. the terms 
breaking continuous translational symmetry of Xri then 
one may write Xr = ko • r/2n + Sx r , with Sxr describing 
fluctuations around the putative average tilt ko deter- 
mined as in Sec. lIIII bv minimizing with respect to ko the 
non-monopole terms in the action. One should regard 
ko as the continuously varying "tilt" the system would 
have were there no cosine (monopole) terms. We will see 
that the true tilt, i.e. (A%) = k, is close but not gen- 
erally equal to ko. The fluctuations of 8xr will then be 
described by a theory of the form 

Suit = [drJ2 U\drSxr\ 2 + p\AS Xr \ 2 } 

— A w cos[27rn.5xr + Tfrkp • r], (33) 

n 

Note that, unlike at the critical point, the fluctuations of 
Sxr at the quadratic level have a non-vanishing (renor- 
malized) stiffness p. Hence the fluctuations of Sxr will be 
bounded, and any non-oscillatory cosine term breaking 
the continuous translational symmetry of Sxr will "pin" 
it, however weak. This pinning corresponds to confine- 
ment in the original gauge theory, and a VBS phase in 
the dimer model. We note in passing that, actually, de- 
pending upon the value of ko and anisotropics in the 
original action, p can be replaced by a more general ten- 
sor. Again, this complication does not modify any of the 
qualitative results of this section, and so will be ignored. 
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For a "generic" value of ko all the cosines oscillate since 
k ■ r will be an irrational multiple of 2ir. There are, how- 
ever, an infinite dense set of values of k for which k • r is a 
rational multiple of 2ir (for all r). In this case, there will 
be some minimal value of n for which the A n term does 
not oscillate (for 8\r — 0). Clearly, at these special val- 
ues of k, this cosine term is relevant and the system is in 
some confined VBS phase. Furthermore, any irrational 
ko is arbitrarily close to one of these rational k values, 
so that, although the cosines in general oscillate on the 
lattice, some of them oscillate extremely slowly. Since in 
general the cosines are not infinitesimally weak (i.e. the 
A„ are finite and non-zero), a sufficiently long wavelength 
oscillation of the cosine term could potentially "pin" the 
5xr field even in such cases. To see whether this occurs, 
let us suppose the n th cosine term oscillates weakly, i.e. 
e mko-r = e *n5k-r for all lattice vec t rs r, with n\Sk\ < 2tt. 
Then, keeping only this cosine term, we have approxi- 
mately 

Stat = [drJ2 U\drS Xr \ 2 +p\AS Xr \ 2 } 

— A n cos[2irn5xr + n 3k ■ r ]- (34) 

Since the fluctuations of Sx r are bounded, we may es- 
timate the effects of the cosine term by ignoring these 
fluctuations and minimizing the action. The minimal 
action configurations are clearly constant in imaginary 
time, d T Sxr — 0. Roughly, then, the field Sx r can either 
be constant, minimizing the stiffness term but gaining no 
(lowering of the) energy from the cosine, or it can choose 
to tilt slightly to take advantage of the cosine term, cost- 
ing some energy from the stiffness. Dimensionally, the 
latter will be favorable if A„ > p|<5k| 2 . Now \Sk\ can be 
made arbitrarily small by increasing n (typically decreas- 
ing as 1/n) , so this inequality will be satisfied if the A„ 
do not decrease too rapidly (i.e. faster than 1/n 2 ) with 
n. However, it is perfectly conceivable that the A n terms 
do decrease faster than 1/n 2 , and in this circumstance, 
there will be incommensurate values of ko for which the 
X field remains unpinned, and all monopole terms re- 
main irrelevant. This has been argued to be the case 
in the immediate neighborhood of RK points in Ref. 0. 
Note that even so, there are commensurate pinned states 
arbitrarily close to this incommensurate state. The self- 
similar succession of various commensurate and incom- 
mensurate states is the incomplete Devil's staircase men- 
tioned above. The term "incomplete" indicates that the 
incommensurate unpinned regions exist (and have finite 
measure, as can also be argued). 

The same arguments apply to the neighborhood of the 
VBS transitions discussed here, and we will sketch the 
reasoning in order to make a few more observations. Note 
that in these cases (e.g. for the honeycomb bilayer), the 
Lifshitz point itself is generic, i.e. can be potentially 
observed in a physical system by varying only one pa- 
rameter. 

First, we comment on a minor subtlety. A naive anal- 



ysis of the continuum field theory, Eq|3 would suggest 
that on the tilted side of the Lifshitz point, the tilt in- 
creases smoothly from zero. This appears to be the case 
since at the Lifshitz point, all cosine (monopole) opera- 
tors are irrelevant. One the tilted side of the transition, 
in fact the tilt does not increase smoothly, but in the 
staircase fashion. This occurs because irrelevant cosine 
operators at the Lifshitz point become relevant on the 
tilted side of the "transition", i.e. these operators are 
dangerously irrelevant in renormalization group parlance. 
In fact, the naive smooth increase in slope (background 
electric field of the gauge theory) occurring on this side 
of the transition is replaced by a slope which contains 
piecewise constant and incommensurate regions, which 
become more and more closely spaced as the Lifshitz 
point is approached, forming an infinite sequence that 
approximates the naive continuous curve (of e.g. |k| ver- 
sus p) arbitrary well if one looks arbitrarily close to the 
Lifshitz point. 

A full description of the devil's staircase is beyond the 
scope of this paper. It is instructive and indicative of the 
general structure to consider some simple "families" of 
plateaus in the tilt k that obtain near the critical point. 
In general, the condition for a plateau is that cos(nk • r) 
does not oscillate on the dual hexagonal lattice. This 
condition is equivalent to requiring that nk is a recipro- 
cal lattice vector of the hexagonal lattice. An arbitrary 
reciprocal lattice vector can be written as nibi + n 2 b 2 , 
with bi = (2tt, 2tt/\/3), b 2 = (0, -4tt/V3). Hence, the 
condition on k is 

k = (£) bl + (^) b 2 (35) 

where n, m, n 2 are integers, and the pairs ni, n and n 2 , n 
can be taken to be relatively prime. The "strongest" 
such commensurate tilts are those with minimal n, since 
these correspond to n-monopole events, which become 
less relevant as n increases. 

In understanding the behavior for small ko , we need to 
investigate those commensurate tilts for which k is small 
but non-zero. Clearly for Eq. [351 this will occur for large 
n. Since the two vectors bi, b 2 are linearly independent, 
the coefficients of both must be small for k to be small. 
Since the two numerators in these coefficients are inte- 
gers, for a given small magnitude |fc| <C 1, one clearly 
then needs at least n > ^S> 1. Larger values of n 

can also yield the same |A;|, by increasing the numerators 
accordingly. However, the largest plateaus in tilt (com- 
ing from the most relevant cosines with minimal n for a 
given |fc|) will be those corresponding to n ~ 

Systems with ko sufficiently near each of these values 
will be pinned and form a "plateau" in k. How wide is 
this plateau? Let us suppose the putative tilt ko ~ k. 
By scaling, the correlation length £ ~ l/l&ol ~ n 3* 1 In 
this situation, we must account for the renormalization of 
the relevant cosine term by the fluctuations on scales less 
than £. From standard renormalization group methods, 
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one expects the renormalized coefficient 

An ~ A„r" 2A , (36) 



where A = -k / (2y K) is the scaling dimension of the one- 
monopole term. Note that, from this reasoning, the A n 
decrease extremely rapidly with n, hence from the above 
argument, incommensurate phases are possible. 

In addition, for large £, the renormalized stiffness is 
small, p ~ p/S, 2 - On scales longer than £, the x field 
is essentially non-fluctuating, so further renormalization 
can be neglected. The criteria for the system to be pinned 
at k is then A„ > p|5fc| 2 . Hence the width of the plateau 



(37) 



Clearly these commensurate plateaus are very narrow 
near the Lifshitz point. Similar estimates were obtained 
in Ref. @. 

Each of the commensurate VBS phases breaks the dis- 
crete translational symmetry of the lattice, and thus must 
undergo a symmetry-restoring transition as temperature 
is increased. At non-zero temperature T, these states will 
be truncated to those commensurate VBS phases whose 
T c s are larger than T. In particular, consider a com- 
mensurately tilted VBS state driven by the n-monopole 
fugacity in the "center" of its plateau, i.e. for k = ko- In 
this case, we may take (5k = in Eq. I|34|l . and the sys- 
tem is described simply by a commensurate sine- Gordon 
model. At T > 0, we may neglect all but the zero Mat- 
subara frequency mode and set d T \ = 0, encompassing 
the quantum effects by using the renormalized parame- 
ters p, A„ above. The symmetry-restoring transition is 
simply the roughening transition of this classical two di- 
mensional sine- Gordon model, which has 



kuT r 



2p 



n 2 e 



(38) 



For the strongest plateaus, recall £ ~ 1/|&| and n ~ 1/|A;|, 
hence one has for these plateaus fcsT c ~ |fc| 4 . The other 
plateaus have even smaller critical temperatures. Since 
all the critical temperatures vanish rapidly as |fc| — ► 0, the 
infinite set of plateaus in the devil's staircase is replaced 
by a finite subset at any non-zero temperature. 

One may also estimate the excitation gap in the associ- 
ated commensurate VBS phase, simply by expanding the 
sine-Gordon term to produce a (Sx) 2 "mass" term. This 

gives a gap E^ BS ~ nV%i ~ ny/X^^ n2A ^ 2 . It may be 
somewhat surprising that a state with an exponentially 
small gap can have such a relatively high critical tem- 
perature (power law in 1/n from Eq. (|38|l ). The physics 
of this is that the excitation with an exponentially small 
gap consists of only small fluctuations of Sx (and hence 
the gauge electric field) which do not perturb the long- 
range order of the VBS state. Indeed, since the VBS 



phases are states of discrete broken symmetry, the ex- 
citations which do disturb this order by connecting the 
different symmetry-related ground states are "droplet" - 
like. In the sine-Gordon language, such a droplet may 
be thought of as a domain wall in Sx, which is wrapped 
around to form a compact "island" inside which Sx is 
shifted by ±27r/n. In particular, the proliferation of such 
thermally excited droplets ultimately will "depin" the 5x 
field and destroy the VBS order above some temperature. 
The minimal radius of such an island is the domain wall 
width itself (since this is larger than the other natural 
cutoff, the correlation length). The energy for such a 
droplet is therefore given by integrating the exponentially 
small sine-Gordon term term over an exponentially large 
area of the size of the domain width (in reality there is 
also a comparable contribution from the p term). These 
two factors compensate to give the relatively large energy 
determining T c . The upshot of this argument is that the 
states with exponentially small excitation gap near the 
zero tilt QCP are not associated with the typical clas- 
sical droplet excitations of a VBS state, but rather are 
evidence of the gapless photon mode obtaining precisely 
at the QCP. 

A comment on the above discussion is in order. Within 
the sine-Gordon treatment, at any temperature above the 
"roughening" temperature, Sx behaves as a free scalar 
field, and vertex operators exp 2ninSx exhibit power law 
correlations. Ultimately, this can be tracked down as 
an artifact of the pure gauge theory. In particular, any 
matter fields included in the model, even gapped ones, 
correspond as discussed earlier to vortices in the x ( or 
Sx) field. At sufficiently high temperatures, these vor- 
tices will certainly unbind. However, for n > 4, it is 
known that such sine-Gordon theories exhibit a "float- 
ing" phase in which power law correlations persist, above 
the roughening temperature T c and below the Kosterlitz- 
Thouless temperature Tkt above which vortices unbind. 
Here UbTkt = p/(8n). Note that one has then 



Tki 



forra > 4. 



(39) 



Hence the VBS phases very close to the Lifshitz point will 
have long-range VBS order only at very low temperatures 
T < T c ~ | A: | 4 but quasi-long-range VBS order up to 
much higher temperatures T < Tkt ~ |fc| 2 - Moreover, 
in the region with quasi-long-range order, there are no 
plateaus in the tilt. 

In gauge theory language the roughening transition 
at T c may be associated with the thermally driven de- 
confinement transition of pure gauge theories. Indeed, 
electric field correlators in this rough phase (gradients of 
the height field) fall off as the inverse square of the dis- 
tance - so this phase may be thought of as a 'thermal 
Coulomb phase'. The gauge charged spinons are loga- 
rithmically interacting in this phase and are bound into 
gauge neutral pairs until they ionize at Tkt leading to a 
plasma of gauge charge that destroys the long range elec- 
tric field correlations of the thermal Coulomb (or rough) 
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phase. This transition is also studied in detail in£, in the- 
ories with a noncompact gauge field and SU(2) symmetric 
spion fields. While the thermal deconfinement transition 
of gauge theories is generally expected only in the ab- 
sence of matter with unit gauge charge (spinons), in two 
spatial dimensions the logarithmic form of the Coulomb 
interactions is strong enough to bind the spinons into 
gauge neutral pairs and hence the transition survives the 
inclusion of spinons. 



IV. THE SINGLE LAYER HONEYCOMB 
LATTICE 

In contrast to the situation analysed above for the bi- 
layer honeycomb aniferromagnet, in a single layer the 
appropriate lattice sine Gordon model has non-zero off- 
sets a r for the x fields on the three sublattices of the 
dual triangular lattice. Thus as explained in Ref. the 
Xr — * ~Xr transformation becomes a symmetry only 
when combined with inversion or a 7r/3 rotation. This 
then leads to the possibility of a cubic invariant in the 
sine Gordon action which drives the transition first order. 

It is interesting to ask about the situation with spin- 
3/2 antiferromagnets on a single layer honeycomb lattice. 
In this case there are no offsets for the Xr in the sine Gor- 
don description of the paramagnetic phase. Nevertheless 
Xr — * ~Xr symmetry (without inversion or 7r/3 rotation) 
is not expected as an exact symmetry of the action. This 
is not required by any of the microscopic symmetries of 
the underlying lattice spin model. Hence we expect that 
a cubic term will still be allowed in the field theory and 
a first order transition will result. 



V. THE SQUARE LATTICE 

We now perform the same analysis for the transition 
in the case of the spin one half quantum antiferromagnet 
on the square lattice. In contrast to the situation on 
the honeycomb lattice, we will find no generic continuous 
transition between VBS states (in both single and bilayer 
cases). 

As before we consider a generic continuum theory that 
is consistent with all the lattice and internal symmetries. 
For the square lattice this takes the form 

S = S c + Si n t + S, non (40) 

Sc = \ J{d TX ? + K[\V\f + Aa{dlx){dlx)W) 



Si, 



'W*xf + {d v xf} + 7;(dxX) 2 {d v xfm 



with a > — 1. The isotropic point corresponds to a = 
and u = v. The last term (S mon ) represents quadru- 
pled monopole events. Note in particular the presence of 
the couplings <r, v which are allowed by the square lat- 
tice symmetry. Consider the critical theory given by the 



quadratic piece of the action S c and determine its sta- 
bility against the inclusion to small monopole tunneling. 
In this case we need to consider quadrupled monopoles, 
as discussed in section [H] and once again it is possible 
to find a range of K, a for which monopole tunneling is 
irrelevant. 

Next, one must consider the stability of the Gaussian 
fixed point described by S c to turning on the quartic 
interaction l|42|) . Again, we perform a one loop renor- 
malization group analysis to determine the fate of these 
couplings. The result of integrating the high wavevec- 
tor modes A(l — dl) < \k\ < A and rescaling, are the 
following RG equations: 



dU 

~~dl 
dV 

IT 



-3[9aU 2 + 2UV + aV 2 }, 
-9[U 2 + 2aUV + V\ 



(43) 
(44) 



where we have used the scaled variables U — uA^ 1 , V — 
vA^ 1 , and the scale factor A = (1 + a + \/l + (t)32itK2 . 
At the isotropic point, these equations are identical to 
(|32|l . and preserve u — v. The anisotropy of the quadratic 
part of the action is present in the parameter a, which is 
unity at the isotropic point but otherwise is given by: 



1 - y/T 



_ 1 2a 
a ~ 3 [_ 'sJTT^-l 



(45) 



thus a £ [1/3, oo). 

We now analyze the RG equations I|43I44|) . and show 
that they imply that the critical Gaussian theory S c is 
generally unstable in the presence of the two quartic op- 
erators. If these were stable critical points, then there 
should be a region in the u, v plane where the flows end 
up at the origin. However we will show that there is 
only a single line in the entire it, v plane, where the cou- 
plings flow into the origin. This implies that stability is 
only attained on a set of measure zero points. For a gen- 
eral choice of the quartic coupling, the flows run away to 
large negative values of u suggesting a first order transi- 
tion. In order to show this property of the RG equations 
(|43I44|I . we construct a function E(U, V) that is invariant 
along any RG trajectory i.e. ^M- = 0. Contour lines of 
this function then represent the RG flows, and we will see 
that there is only a single contour connected to the origin. 
This function is most conveniently represented in terms 
of the rotated coordinates, U+ = all + V, U- = U — aV. 
Then, 



E 



{a 2 + 3)U' 2 + AaU+U- + (3a 2 + l)Ul ' 



(46) 



which can be checked to be invariant under the RG flows. 
It may be seen that the origin must correspond to E = 
0, by approaching it in any direction. The only other 
points for which the invariant function vanishes is the 
line U- = 0, which corresponds to u — av. These are 
the only points that could flow to the origin. Apart from 
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this set of measure zero, none of the other points in the 
u, v plane ever reach the origin under the RG flows but 
rather flow to regions with large negative values of u. 
Thus, for the case with square symmetry, the Gaussian 
critical theory is generically unstable - the flows suggest a 
first order transition in the absence of special fine-tuning. 



VI. RK POINTS 

The quantum dimer model Hamiltonians studied 
were shown to have a special point - the RK point - at 
which the wavefunction is an equal superposition of all 
dimer configurations. Equal time correlation functions 
can then be evaluated from the classical dimer model, 
which has been extensively studied^. In this section we 
address how those results fit into the framework discussed 
here - for the case of the single layer bipartite lattices. 
For example we may ask in the case of the honeycomb 
lattice where a transition through a multicritical point 
can obtain by tuning two parameters, whether the RK 
wavefunction corresponds to the the ground state wave 
function of any of these fixed points. In fact, we will con- 
clude that while the RK points in both the square and the 
honeycomb lattice cases can reproduce critical properties 
of some point on the line of fixed points obtained after 
tuning a few parameters in the generic models, they rep- 
resent very special multicritical points in terms of their 
position in the phase diagram of these generic models. 
Thus accessing these points requires fine tuning an infi- 
nite number of independent operators, which do not af- 
fect the critical properties but change the phase structure 
in the immediate vicinity of the point (dangerously irrel- 
evant operators). This fine tuned nature of the RK point 
can immediately be seen by noting that at the RK point, 
the ground state wavefunction in each winding number 
(tilt) sector has exactly the same energy. Reproducing 
this democratic treatment of all winding number sectors 
within a height model representation of the RK point will 
require tuning an infinite number of parameters to zero 
in the bare Hamiltonian - even though they may be as- 
sociated with operators that are (dangerously) irrelevant 
at the critical point. 

If the RK point requires fine tuning infinitely many 
parameters, one may ask how it is accessed so readily 
in the quantum dimer model. The reason is that the 
dimer model usually contains only single plaquette terms, 
where the RK point is accessed by tuning just one param- 
eter. However, including processes that involve several 
plaquettes will require tuning increasingly larger number 
of parameters to obtain the equal dimer superposition 
of the RK point. Generically the dimer representation 
of spin systems will include processes that involve arbi- 
trarily large number of plaquettes. Although these decay 
with increasing size, they are generally non vanishing, 
and infinite fine tuning will be then be required to reach 
the RK point. 

For the square lattice, the field theory that reproduces 





Tilt=0 




\ C/P 
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Staggered 




(c) 



FIG. 2: Schematic depiction of the phase diagram of VBSs, 
the vertical axis in the plots is roughly the the parameter p 
and the staggered state has the maximum tilt. The generic 
phase diagrams expected for the bilayer honeycomb lattice 
are shown in (a) and (b) (and also of the single layer honey- 
comb lattice after tuning one parameter, the cubic term, to 
zero), (a) The continuous transition (shown with the dashed 
line) is from a zero tilt phase to a region where the tilt (Q) 
exhibits a a devil's staircase structure, contours of equal tilt 
(thin solid lines) are shown. The critical line ends in a mul- 
ticritical point M' beyond which the transition is first order 
(solid line). The horizontal axis here may be thought of as 
the coefficient of the quartic term, (b) An alternate scenario, 
there is again a continuous transition to a region with the 
devil's staircase structure for the tilt. Here however, there 
is a multicritical point M that is adjacent to the staggered 
state which could control a zero tilt to staggered state tran- 
sition. The horizontal axis here may be thought of as the 
energy cost of the maximally tilted (staggered) state, (c) The 
RK point for the single layer square/honeycomb lattice - an 
infinite number of parameters need to be tuned to access this 
plane. Exact degeneracy of the different winding number sec- 
tor ground states implies that states with arbitrary tilt lie 
infinitesimally close to the RK point as shown. 



asymptotic properties of the RK point was studied by 
Henley^. It was found there that a quadratic action l|40|) 
with K = 7r 2 /4, (7 = gives the long distance dimer cor- 
relations at the RK point (amusingly this is the parame- 
ter value at which the two monopole insertion operator is 
marginal) . The fact that the ground states in all winding 
number sectors are degenerate at the RK point implies 
that it is possible to access states with arbitrary values 
of the tilt my moving infinitesimally away from the RK 
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point. As a result, the topology of the phase diagram 
near the RK point is as shown in figure J2t) . An infinite 
number of parameters have been set to zero to access this 
plane of the phase diagram. While the continuous transi- 
tion from the zero tilt state leads to the devil's staircase 
of tilted states, a direct transition from a zero tilt state 
to the staggered state can be made by crossing the RK 
point. Thus the RK point terminates the critical line, 
and sits on a line of first order transitions. The special 
nature of the RK point, even as a multicritical point, 
can be seen by comparing it to the generic multicritical 
points of the bilayer honeycomb lattice (which require 
tuning of two parameters to reach) denoted M', M in 
figures [2t,b. (These can also be thought of as higher 
order multicritical points in the single layer honeycomb 
model, which will require fine tuning of an additional pa- 
rameter to reach the plane depicted in the diagrams.) In 
the first scenario depicted in figure Ek, the multicritical 
point M' terminates both the first order transition line 
and the critical line, but it does not allow for a direct 
transition from the zero tilt state to the staggered state. 
In the other scenario figure JSb) the generic multicritical 
point M does allow for a direct transition from a zero tilt 
state to the staggered state. However, it is fundamentally 
distinct from the RK point, as can be seend from the dif- 
ference in topology of the phase structure around the M 
and RK points, which may be characterized as follows. 
We ask what phases may be accessed from these points 
by a small change of bare parameters. For the generic 
multicritical point M, the staggered state, the zero tilt 
state, and a state with infinitesimal tilt can be accessed. 
For the RK point however, besides the staggered and 
zero tilt state, states with arbitrary values of the tilt can 
also be accessed in this manner, as shown in the figure. 
This follows from the exact degeneracy of ground states 
in different winding number sectors at the RK point. 



VII. CONCLUSIONS 

In this paper we have shown that at least in certain 
instances there are direct second order transitions be- 
tween distinct valence bond solid phases. Deconfincment 
obtains at the critical point though both phases are con- 
ventional and confined. More precisely the critical theory 
may be viewed as a gapless U(X) gauge theory with irrel- 
evant instantons. The spin carrying excitations in either 
phase are gapped spin-1 magnons or their composites. 
Right at the critical point however the spin gap does 
not close but the magnons decay into (gapped) spin-1/2 
spinons. 

All of this structure is very similar to the other exam- 
ples of 'deconfined' quantum criticality studied in Ref. 2 . 
However there are some interesting differences in the de- 
tails. First the critical points discussed in the present 
paper are controlled by a critical fixed line with contin- 
uously varying exponents. Second (upto a marginally ir- 
relevant non-linear operator) all points on this fixed line 



have a simple free field description. The emergence of a 
(topological) global U(l) symmetry that seems generic to 
deconfined quantum criticality obtains for the transitions 
in this paper as well. However the free field description 
implies an enormous number of further emergent symme- 
tries which are specific to these transitions. Third, the 
dynamic scaling exponent z is 2 in the present example 
(compared to z = 1 at the Neel-VBS transition). Finally 
as detailed in Section Till CI there is rich and interesting 
structure with an infinite number of transitions on one 
side of the deconfined critical point. 

One off-shoot of these results is a clarification of the 
place of the solvable RK point of quantum dimer models 
in a more general context of phase transitions in quantum 
magnets. We find that the RK point corresponds to a 
special multicritical point. 

We also showed that these interesting phase transitions 
are best realized in bilayer spin-1/2 honeycomb lattice 
quantum antiferromagnets. It would be interesting for 
numerical work to explore specific spin models on such 
bilayers where these transitions can be accessed. 

Several extensions of our results are possible. It should 
be possible to examine the role of various perturbing 
fields at the critical point as well as the effects of finite 
temperature. It should also be readily possible to exam- 
ine transitions between different VBS phases in higher 
spin quantum magnets. We leave these for future work. 

Since the original submission of an electronic preprint 
of this work, Ref. |^ appeared which considered some of 
the same questions especially those regarding the single 
layered systems discussed here. They correctly pointed 
out the first order nature of the transition on the single 
layer honeycomb lattice, and the 'incompleteness' of the 
devil's staircase - which were the points of disagreement 
with the earlier version of this work. These points have 
been corrected and briefly mentioned in the present work. 
However in contrast to Ref. we have focused here 
on the case of the bilayer honycomb quantum magnet 
where a generically continuous transition between VBS 
states is realised, and have studied some of the interesting 
properties of this generic quantum critical point. 
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APPENDIX A: CRITICAL GAUGE THEORY 

In this Appendix we will explicitly display the form 
of the continuum action of the non-compact U(X) gauge 
theory that corresponds to the critical sine Gordon the- 
ory. A general continuum Hamiltonian for a non-compact 
U(l) gauge theory in D = 2+1 dimensions takes the form 

H = J d 2 xp\E\ 2 + K\V x E| 2 + B 2 + ■ ■ ■ (Al) 

together with the Gauss law constraint 

V • E = 0. (A2) 

Here E refers to the electric field while B = Z'VxA = 
eijdiAj is the magnetic field. We work in the Coulomb 
gauge V- A = 0. As usual the components of the electric 
field Ei (i = x, y) and the (transverse) components aj of 
the vector potential are canonically conjugate: 

[^(x), A,-(x')] = -%P tj (x - x'), (A3) 

k k 

where Pij has Fourier components Sij — 1 and projects 
out the transverse component. 

The Gauss law constraint may be solved by writing 

Ei = Cijdjx, (A4) 



where X 1S a scalar field. Assume that the commutator 
of x an d Aj takes the form 

[ X (x),A i (x / )]=t/ i (x-x'). (A5) 

The correct commutator between E and A is reproduced 
if we impose 

dkdkfj (x - x') = -Py (x - x'). (A6) 

This then implies the commutator 

[ X (x),B(x')]=^(x-xO, (A7) 

so that the magnetic field B is conjugate to X - The 
Hamiltonian may now be rewritten 

H = J d 2 xp\V X \ 2 +K(X7 2 x ) 2 +B 2 (A8) 

As B and x are canonically conjugate we reproduce 
the continuum free field action of the sine Gordon the- 
ory. Remembering that the critical theory has p = 
we can immediately read off from Eqn. IA1I the contin- 
uum gauge theory Hamiltonian that describes the criti- 
cal point. This is readily diagonalized explicitly to find 
a quadratic dispersing photon. 

Amusingly, at p = 0, the Hamiltonian (|A8|) exhibits 
a kind of self 'duality', obtained by exchanging the roles 
of B and V 2 %- That is, if we introduce the field <j> such 
that B = ^t=^, and its conjugate field = — \/jfV 2 x, 
then one obtains the same critical Hamiltonian (|A8l with 
p = 0) but with K^l/K and (B, x) -> (n , <f>). 
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